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Q ! Quarkonia play an important role in several high energy experiments. The diversity, 
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quantity and accuracy of the data still under analysis and currently being collected is 
impressive and includes: data on quarkonium formation from BES at the Beijing Elec- 



ts ■ tron Positron Collider (BEPC), E835 at Fermilab, and CLEO at the Cornell Electron 

^ . Storage Ring (CESR); clean samples of charmonia produced in B-decays, in photon- 

', photon fusion and in initial state radiation, at the B-meson factories, BaBar at PEP-II 

I I and Belle at KEKB, including the unexpected observation of large amounts of asso- 

dated (cc) (cc) production and the observation of new and possibly exotics quarkonia 
1^ • states; the CDF and DO experiments at Fermilab measuring heavy quarkonia produc- 

tion from gluon-gluon fusion in pp annihilations at 2 TeV; the Selex experiment at 

* I 

^ , Fermilab with the preliminary observation of possible candidates for doubly charmed 

H I baryons; ZEUS and HI, at DESY, studying charmonia production in photon-gluon 

fusion; PHENIX and STAR, at RHIC, and NA60, at CERN, studying charmonia 
production, and suppression, in heavy-ion collisions. This has led to the discovery of 
new states, new production mechanisms, new decays and transitions, and in general 
to the collection of high statistics and precision data sample. In the near future, even 
larger data samples are expected from the BES-III upgraded experiment, while the 
B factories and the Fermilab Tevatron will continue to supply valuable data for few 
years. Later on, new experiments at new facilities will become operational (the LHC 
experiments at CERN, Panda at GSI, hopefully a Super-B factory, a Linear Collider, 
etc.) offering fantastic challenges and opportunities in this field. A comprehensive 
review of the experimental and theoretical status of heavy quarkonium physics may 
be found in the Cern Yellow Report prepared by the Quarkonium Working Group 
[1]. Many excellent reviews of the field have been presented at this meeting 
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On the theory side, systems made by two heavy quarks are a rather unique lab- 
oratory. They are characterized by the existence of a hierarchy of energy scales in 
correspondence of which one can construct a hierarchy of nonrelativistic effective 
field theries (NR EFT), each one with less degrees of freedom left dynamical and thus 
simpler. Some of these physical scales are large and may be treated in perturbation 
theory. The occurrence of these two facts makes two heavy quark systems accessible 
in QCD. In particular, the factorization of high and low energy scales realized in 
the EFTs allows us to study low energy QCD effects in a systematic way. Today 
the remarkable progress in the construction of these nonrelativistic EFTs together 
with the advance in lattice QCD give us well based theory tools to investigate heavy 
quarkonia. 

Therefore, on the one hand the progress in our understanding of EFTs makes it 
possible to move beyond phenomenological models and to provide a systematic de- 
scription from QCD of several aspects of heavy-quarkonium physics. On the other 
hand, the recent progress in the measurement of several heavy-quarkonium observ- 
ables makes it meaningful to address the problem of their precise theoretical deter- 
mination. In this situation heavy quarkonium becomes a very special and relevant 
system to advance our understanding of strong interactions and our control of some 
parameters of the Standard Model. 

Here I will briefly review some of the recent developments in the construction of 
NR EFTs with the main emphasis on the physical applications. For some reviews see 

lailElE]. 

2 Scales and Effective Field Theories 

The description of hadrons containing two heavy quarks is a rather challenging prob- 
lem, which adds to the complications of the bound state in field theory those coming 
from the nonperturbative QCD low-energy dynamics. A simplification is provided 
by the nonrelativistic nature suggested by the large mass of the heavy quarks and 
manifest in the spectrum pattern. As nonrelativistic systems, quarkonia are charac- 
terized by three energy scales, hierarchically ordered by the heavy quark velocity in 
the center of mass frame f <C 1: the mass m (hard scale), the momentum transfer mv 
(soft scale), which is proportional to the inverse of the typical size of the system r, 
and the binding energy mv'^ (ultrasoft scale), which is proportional to the inverse of 
the typical time of the system. In bottomonium f ^ ~ 0.1, in charmonium v"^ ~ 0.3, in 
ti V ^ 0.15. In perturbation theory v ~ Og. Feynman diagrams will get contributions 
from all momentum regions associated with these scales. Since these momentum re- 
gions depend on as, each Feynman diagram contributes to a given observable with 
a series in ag and a non trivial counting. Besides, the ctg associated to different mo- 
mentum region are evaluated at different scales. For energy scales close to Aqcd, the 
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scale at which nonperturbative effects become dominant, perturbation theory brealcs 
down and one has to rely on nonperturbative methods. Regardless of this, the non- 
relativistic hierarchy m ^ mv ^ mv'^, m ^ Aqcd will persist also below the Aqcd 
threshold. 

The wide span of involved energy scales makes also a lattice calculation in full QCD 
extremely challenging. However, it is possible to exploit the existence of a hierarchy 
of scales by introducing a hierarchy of nonrelativistic effective field theories. Lower 
energy EFTs may be constructed by systematically integrating out modes associated 
to energy scales not relevant for the two quark system. Such integration is made 
in a matching procedure that enforces the equivalence between QCD and the EFT 
at any given order of the expansion in v. Any prediction of the EFT is therefore a 
prediction of QCD with an error of the size of the neglected order in v. By integrating 
out the hard modes, one obtains Nonrelativistic QCD P [TOl [13]. In such EFT, soft 
and ultrasoft scales are left dynamical and still their mixing complicates calculations 
and power counting. In the last few years the problem of systematically treating 
the remaining dynamical scales in an EFT framework has been addressed by several 
groups |Tl] and has now reached a good level of understanding. So one can go down 
one step further and integrate out also the soft scale in a matching procedure to 
the lowest energy and simplest EFT that can be introduced for quarkonia, where 
only ultrasoft degrees of freedom remain dynamical. Here I will review potential 
NRQCD [71 E] , for an alternative and equivalent EFT (in the case in which Aqcd is 
the smallest scale) see [12]. In the case in which the soft scale is of the same order of 
Aqcd, the matching to pNRQCD is still possible but it is nonperturbative. 

3 NonRelativistic QCD (NRQCD) 

NRQCD [9} [To] is the EFT for two heavy quarks that follows from QCD by integrating 
out the hard scale m. Only the upper (lower) components of the Dirac fields remain 
relevant for dymanical quarks (antiquarks) at energies lower than m. Thus quark and 
antiquarks are described in terms of two-components Pauli spinor fields. The part 
of the NRQCD Lagrangian bilinear in the heavy quark fields is the same as Heavy 
Quark Effective Field Theory (HQET) (for a review see [I4j) but for the case of two 
heavy quarks also four fermion operators have to be considered. The Lagrangian is 
organized as an expansion in v and as{m): 

>Cnrqcd = '^Cn{m,fi) X 0„(/i,mt;,mt;^ Aqcd)M"- (1) 

n 

The NRQCD matching coefficients c„ are series in ag and encode the ultraviolet 
physics that has been integrated out from QCD. The low energy operators 0„ are 
constructed out of two or four heavy quark/antiquark fields plus gluons. They are 
counted in powers of v. Since two scales, soft and the ultrasoft, are dynamical, the 
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power counting in v is not unambiguous. The imaginary part of the coefficients of 
the 4-fermion operators contains the information on heavy quarkonium annihilations. 
The NRQCD heavy quarkonium Fock state is given by a series of terms, increasingly 
subleading, where the leading term is a QQ in a color singlet state and the ffist 
correction, suppressed in v, comes from a QQ in a color octet state plus a gluon. The 
NRQCD Lagrangian can be used for studies of spectroscopy (on the lattice), inclusive 
decays and electromagnetics threshold production of heavy quarkonia. 

4 potential NonRelativistic QCD (pNRQCD) 

pNRQCD [3 El [3] is the EFT for two heavy quark systems that follows from NRQCD 
by integrating out the soft scale mv. Here the role of the potentials and the quan- 
tum mechanical nature of the problem are realized in the fact that the Schrodinger 
equation appears as zero order problem for the two quark states. We may distinguish 
two situations: 1) weakly coupled pNRQCD when mv ^ Aqcd, where the matching 
from NRQCD to pNRQCD may be performed in perturbation theory; 2) strongly 
coupled pNRQCD when mv ~ Aqcd, where the matching has to be nonperturbative. 
Recalling that r"^ ~ mv, these two situations correspond to systems with inverse 
typical radius smaller than or of the same order as Aqcd- 

4.1 Weakly coupled pNRQCD 

The effective degrees of freedom that remain dynamical are: low energy QQ (Pauli 
spinor) states that can be decomposed into a singlet field S and an octet field O 
under colour transformations, have energy of order Aqcd, ttt-v'^ and momentum p of 
order mv; low energy (ultrasoft (US)) gluons A^(R, t) with energy and momentum 
of order Aqcd, mv"^. All the gluon fields are multipole expanded (i.e. expanded in 
the quark-antiquark distance r, R being the center of mass). The Lagrangian is then 
given by terms of the type 

^^^^ X K(r/i', r/.) X 0„(/i', mv\ Aqcd) r". (2) 
m" 

where the matching coefficients are inherited from NRQCD and contain the logs in 
the quark masses, the pNRQCD potential matching coefficients Vn encode the non- 
analytic behaviour in r and the low energy operators 0„ are constructed in terms of 
singlet, octet fields and ultrasoft gluons. At leading order in the multipole expansion, 
the singlet sector of the Lagrangian gives rise to equations of motion of the Schrodinger 
type. Each term in the pNRQCD Lagrangian has a definite power counting. The bulk 
of the interaction is carried by potential-like terms, but non-potential interactions, 
associated with the propagation of low energy degrees of freedom are present as well. 
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Such retardation (or non-potential) effects start at the next-to-leading order (NLO) 
in the multipole expansion and are systematically encoded in the theory and typically 
related to nonperturbative effects [8]. There is a systematic procedure to calcolate 
corrections in v to physical observables: higher order perturbative (bound state) 
calculations in this framework become viable. In particular the EFT can be used 
for a very efficient resummation of large logs (typically logs of the ratio of energy 
and momentum scales) using the renormalization group (RG) adapted to the case 
of correlated scales [12] : Poincare invariance is not lost, but shows up in some 
exact relations among the matching coefficients jl6]. The renormalon subtraction 
may be implemented systematically obtaining a perturbative series better behaved 
and allowing a factorization of the genuine QCD nonperturbative effects. 



In this case the matching to pNRQCD is nonperturbative. Away from threshold 
(precisely when heavy-light meson pair and heavy hybrids develop a mass gap of 
order Aqcd with respect to the energy of the QQ pair), the quarkonium singlet 
field S remains as the only low energy dynamical degree of freedom in the pNRQCD 
Lagrangian (if no ultrasoft pions are considered), which reads [T^ [20t [5]: 



The matching potential Vs{r) is a series in the expansion in the inverse of the quark 
masses: static, 1/m and 1/m? terms have been calculated, see [T9l [20] . They involve 
NRQCD matching coefficients and low energy nonperturbative parts given in terms 
of Wilson loops and field strengths insertions in the Wilson loop. In this regime we 
recover the quark potential singlet model from pNRQCD. However the potentials are 
calculated from QCD in the formal nonperturbative matching procedure. An actual 
evaluation of the low energy part requires lattice evaluation [17] or QCD vacuum 
models calculations [T8l [27]. 

5 Applications 

The condition m <C Aqcd always holds and thus the first matching from QCD to 
NRQCD is a perturbative matching. NRQCD describes in principle all heavy quarko- 
nia states and physical processes. However, since still the soft and the ultrasoft scales 
are dynamical, the power counting is not unambiguous and in some cases may differ 
from the perturbative inspired BBL counting [TUl [Ml [25] ■ The number of nonper- 
turbative operators tend to increase with the order of the expansion in f , and their 
expectation values depend both on the quarkonium states and the US gluons. The 



4.2 Strongly coupled pNRQCD 




(3) 
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NRQCD lattice implementation still requires the calculation of the NRQCD matching 
coefficients in the lattice regularization, which is still missing in many cases. Being 
NRQCD a nonrenormalizable theory at the leading order Lagrangian, NRQCD lattice 
calculations maybe tricky. 

The lowest energy EFT, pNRQCD, is simpler and as such may be more predictive. 
However, in the present formulation, it is valid only for states away from threshold. 
Since we are now integrating out also the soft scale, it is important to establish 
when Aqcd sets in, i.e. when we have to resort to non-perturbative methods. For 
low-lying resonances, it is reasonable, although not proved, to assume mv"^ > Aqcd- 
The system is weakly coupled and we may rely on perturbation theory, for instance, 
to calculate the potential. The theoretical challenge here is performing higher-order 
calculations and the goal is precision physics. For high-lying resonances, we assume 
mv ~ Aqcd- The system is strongly coupled and the potential must be determined 
non-perturbatively, for instance, on the lattice. The theoretical challenge here is pro- 
viding a consistent framework where to perform lattice calculations and the progress 
is measured by the advance in lattice computations. The number of nonperturba- 
tive operators maybe be greatly reduced with respect to NRQCD, since a further 
factorization at the soft scale is realized and nonperturbative contributions become 
typically only a function of the US gluons. The pNRQCD leading order strongly 
coupled Lagrangian is renormalizable allowing in principle a straightforward lattice 
implementation. 

Both in NRQCD and pNRQCD a source of concern may arise from the large 
corrections in the charmonium case. Large (renormalon-like) perturbative contri- 
butions in the matching coefficients need to be properly taken care, resummed and 
subtracted. 

6 QCD potentials 

The QCD potentials achieve a well defined status and definition only in pNRQCD: 
they are the matching coefficients of the EFT and as such there is a well defined 
procedure to calculate them. They depend on the scale of the matching. In weakly 
coupled pNRQCD the soft scale is bigger than Aqcd and so the singlet and octet 
potentials have to be calculated in the perturbative matching. In ^21j a determination 
of the singlet potential at three loops leading log has been obtained inside the EFT 
which gives the way to deal with the well known infrared singularity arising in the 
potential at this order [36] . From this, as in the V regularization can be obtained, 
showing a dependence on the infrared behaviour of the theory at this order and for this 
regularization. The finite terms in the singlet static potential at three loops are not 
yet known but has been estimated [23]. Recently, also the logarithmic contribution at 
four loops has been calculated [23]. The three loop renormalization group improved 
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calculation of the static singlet potential has been compared to the lattice calculation 
and found in good agreement up to about 0.25 fm [22]. The static octet potential is 
known at two loops [25] and again agrees well with the lattice data [26]. 

At a scale fi such that mv ~ Aqcd S> ^ mt>^, confinement sets in and the 
potentials become admixture of perturbative terms, inherited from NRQCD, which 
encode high-energy contributions, and non-perturbative objects. Strongly coupled 
pNRQCD gives us the general form of the potentials obtained in the nonperturbative 
matching to QCD in the form of Wilson loops and Wilson loop chromoelectric and 
chromomagnetic field strengths insertions [T9l [20] , very well suited for lattice calcula- 
tions. These will be in general complex valued functions. The real part controls the 
spectrum and the imaginary part controls the decays. 

The real part of the potential has been one of the first quantities to be calculated 
on the lattice (for a review see [IZj). In the last year, there has been some remark- 
able progress. In [2S], the 1/m potential has been calculated for the first time. The 
existence of this potential was first pointed out in the pNRQCD framework ^19]. A 
1/m potential is typically missing in potential model calculations. The lattice result 
shows that the potential has a 1/r behaviour, which, in the charmonium case, is of 
the same size as the 1/r Coulomb tail of the static potential and, in the bottomo- 
nium one, is about 25%. Therefore, if the 1/m potential has to be considered part of 
the leading-order quarkonium potential together with the static one, as the pNRQCD 
power counting suggests and the lattice seems to show, then the leading-order quarko- 
nium potential would be, somewhat surprisingly, a flavor-dependent function. In [29] , 
spin-dependent potentials have been calculated with unprecedented precision. In the 
long range, they show, for the first time, deviations from the fiux-tube picture of chro- 
moelectric confinement [SHIEZI- The knowledge of the potentials in pNRQCD could 
provide an alternative to the direct determination of the spectrum in NRQCD lattice 
simulations: the quarkonium masses would be determined by solving the Schrodinger 
equation with the lattice potentials. The approach may present some advantages: the 
leading-order pNRQCD Lagrangian, differently from the NRQCD one, is renormaliz- 
able, the potentials are determined once for ever for all quarkonia, and the solution of 
the Schrodinger equation provides also the quarkonium wave functions, which enter 
in many quarkonium observables: decay widths, transitions, production cross- 
sections. The existence of a power counting inside the EFT selects the leading and 
the subleading terms in quantum-mechanical perturbation theory. Moreover, the 
quantum mechanical divergences (typically encountered in perturbative calculations 
involving iterations of the potentials, as in the case of the iterations of spin delta 
potentials) are absorbed by NRQCD matching coefficients. Since a factorization be- 
tween the hard (in the NRQCD matching coefficients) and soft scales (in the Wilson 
loops or nonlocal gluon correlators) is realized and since the low energy objects are 
only glue dependent, confinement investigations, on the lattice and in QCD vacuum 
models become feasible [TS] f27\. 
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The potentials evaluated on the lattice once used in the Schrodinger equation 
produce the spectrum. The calculations involve only QCD parameters (at some scale 
and in some scheme). 



7 Precision determination of Standard Model pa- 
rameters 

Given the advancement in the EFTs formulation and in the lattice calculations as well 
as the existence of several high order perturbative bound state calculations, quarkonia 
may become a very appropriate system for the extraction of precise determination of 
the Standard Model parameters like and the heavy quark masses. Such precise 
determinations are important for physics inside and beyond the Standard Model. 
Inside the QWG (www.qwg.to.infn.it) there is a topical group for such studies and in 
the QWG YR [1\ there is a dedicated chapter. 

7.1 c and b mass extraction 

The lowest heavy quarkonium states are suitable systems to extract a precise de- 
termination of the mass of the heavy quarks b and c. Perturbative determinations 
of the T(IS') and J/ip masses have been used to extract the b and c masses. The 
main uncertainty in these determinations comes from nonperturbative nonpotential 
contributions (local and nonlocal condensates) together with possible effects due to 
subleading renormalons. These determinations are competitive with those coming 
from different systems and different approaches (for the b mass see e.g. [HE]). We 
report some recent determinations in Tab. 17.11 

A recent analysis performed by the QWG [Ij and based on all the previous deter- 
minations indicates that at the moment the mass extraction from heavy quarkonium 
involves an error of about 50 MeV both for the bottom (1% error) and in the charm 
(4% error) mass. It would be very important to be able to further reduce the error 
on the heavy quark masses. 

7.2 Determinations of o^g. 

Heavy quarkonia leptonic and non-leptonic inclusive and radiative decays may provide 
means to extract a^. The present PDG determination of from bottomonium pulls 
down the global as average noticeably |T]. Recently, using the most recent CLEO 
data on radiative T(IS') decays and dealing with the octet contributions within weakly 
coupled pNRQCD, a new determination of as (Mt(i5)) = 0.184t[J:!]i3 has been obtained 
t43j, which corresponds to a^iMz) = O.lldtomt 

in agreement with the central value 

of the PDG |73j and with competitive errors. 
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reference 


order 


fntiiTT-b) (GeV) 


1871 
]P ' J 


NNNLO* 

J. 1 J. 1 J. 1 1 J V, ./ 


4 210 ± 090 ± 025 


[881 


NNLO +charm 


4.190 ±0.020 ±0.025 


I90l 


NNLO 


4 24 ± 10 


[891 


NNNLO* 


4.346 ± 0.070 


m\ 


NNNLO* 


4.20 ±0.04 


m\ 


NNNLO* 


4.241 ± 0.070 


[93J 


NNLL* 


4.19 ±0.06 


reference 


order 


rricijfic) (GeV) 


m 


NNLO 


1.24 ±0.020 


m 


NNLO 


1.19 ±0.11 



Table 1: Different recent determinations of fnh{Tnf,) and Wicirn^ in the MS scheme 
from the bottomonium and the charmonium systems. The displayed results either 
use direct determinations or non-relativistic sum rules. Here and in the text, the * 
indicates that the theoretical input is only partially complete at that order. 

7.3 Top-antitop production near threshold at ILC. 

In |93l [91] the total cross section for top quark pair production close to threshold 
in e±e- annihilation is investigated at NNLL in the weakly coupled EFT. The sum- 
mation of the large logarithms in the ratio of the energy scales is achieved with the 
renormalization group (for correlated scales) and significantly reduces the scale de- 
pendence of the results. Studies like these will make feasible a precise extractions of 
the strong coupling, the top mass and the top width at a future ILC. 

8 Spectra 

The NRQCD Lagrangian is well suited for lattice calculations [21] • The quark prop- 
agators are the nonrelativistic ones and since the heavy-quark mass scale has been 
integrated out, for NRQCD on the lattice, it is sufficient to have a lattice spacing a 
as coarse as m ^ 1/a ^ mv. A price to pay is that, by construction, the contin- 
uum limit cannot be reached. Another one is that the NRQCD Lagrangian has to 
be supplemented by matching coefficients calculated in lattice perturbation theory, 
which encode the contributions from the heavy-mass energy modes that have been 
integrated out. A recent unquenched determination of the bottomonium spectrum 
with staggered sea quarks can be found in [32]. The fact that all matching coefficients 
of NRQCD on the lattice are taken at their tree-level value induces a systematic er- 
ror of order a^v"^ for the radial splittings and of order for the fine and hyperfine 
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splittings. 

Inside pNRQCD we have to consider separately systems with a small interquark 
radius (low-lying states) and systems with a radius comparable or bigger than the 
confinement scale Aq^d (high- lying states). It is difficult to say to which group a spe- 
cific resonance may belong, since there are no direct measurements of the interquark 
radius. Electric dipole transitions or quarkonium dissociation in a medium, once a 
well founded theory treatment of such processes will be given, may give a clear cut 
procedure. At the moment one uses the typical EFT approaches assuming that a 
particular scales hierarchy holds and checking then a posteriori that the prediction 
and the error estimated inside such framework are consistent with the data. 

Low-lying QQ states are assumed to realize the hierarchy: m ^ mv ^ mv'^ > Aqcd 
and they may be described in weakly coupled pNRQCD. 



Be mass (MeV) 


[84J (expt) 




83J (lattice) 


82J (NNLO) [81J (NNLO) [88] (NNLO) 


6287 ±4.8 ± 1.1 


6304 ±121^^ 6326(29) 6324(22) 6307(17) 



Table 2: Different perturbative determinations of the mass compared with the 
experimental value and a recent lattice determination. 

Once the heavy quark masses are known, one may use them to extract other 
quarkonium ground-state observables. The Be mass has been calculated at NNLO in 
[HSIIHIIIHH], see TableO These values agree well with the unquenched NRQCD/Fermilab 
method) lattice determination of [83] , which shows that the Be mass is not very sen- 
sitive to non-perturbative effects. This is confirmed by a recent measurement of the 
Be in the channel Be J/ipir hj the CDF collaboration at the Tevatron; they ob- 
tain with 360 pb"^ of data Mb, = 6285.7 ± 5.3 ± 1.2 MeV [H], while the latest 
available figure based on 1.1 fb~^ of data is Mb, = 6276.5 ± 4.0 ± 2.7 MeV (see 
http: / /www-cdf.fnal.gov/physics/new/bottom/060525.blessed-bc-mass/ ). 

The bottomonium and charmonium ground-state hyperfine splitting has been cal- 
culated at NLL in [95]. Combining it with the measured T(IS') mass, this determi- 
nation provides a quite precise prediction for the r]b mass: M^^ = 9421 ± lOtl MeV, 
where the first error is an estimate of the theoretical uncertainty and the second one 
refiects the uncertainty in as- Note that the discovery of the rif, may provide a very 
competitive source of as at the bottom mass scale with a projected error at the Mz 
scale of about 0.003. Similarly, in [96], the hyperfine splitting of the Be was calculated 
at NLL accuracy: M^. - Mb, = 65 ± 24l}^ MeV. 

High-lying QQ states are assumed to realize the hierarchy: m ^ mv ~ Aqcd ^ 
mf ^. A first question is where the transition from low-lying to high-lying takes place. 
This is not obvious, because we cannot measure directly mv. Therefore, the answer 
can only be indirect and, so far, there is no clear agreement in the literature. A 
weak-coupling treatment for the lowest-lying bottomonium states [n = 1, n = 2 and 
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also for the T{3S)) appears to give positive results for the masses at NNLO in [81] 
and at N'^LO* in [37]. The result is more ambiguous for the fine splittings of the 
bottomonium IP levels in the NLO analysis of [42] and positive only for the T(15') 
state in the N^LO* analysis of [38]. 




Masses of high-lying quarkonia may be accessed either using the lattice nonpertur- 
bative potentials inside a Schrodinger equation [33] or via a direct lattice pNRQCD 
calculation. 

9 Transitions and decays 
9.1 Inclusive Decays 

NRQCD gives a factorization formula for heavy quarkonium inclusive decay widths, 
precisely it factorizes four-fermion matching coefficients and matrix elements of four 
fermion operators [10]. Color singlet operator expectation values may be easily related 
to the square of the quarkonium wave functions (or derivatives of it) at the origin. 
Octet contributions remain as nonperturbative matrix elements of operators evaluated 
over the quarkonium states. In some situations the octet contributions may not be 
suppressed and become as relevant as the singlet contributions in the NRQCD power 
counting. In particular octet contributions may reabsorb the dependence on the 
infrared cut-off of the Wilson coefficients, solving the problem that arised originally 
in the color singlet potential model |62j . 

Systematic improvements are possible, either by calculating higher-order correc- 
tions in the coupling constant or by adding higher-order operators. 

In order to describe electromagnetic and hadronic inclusive decay widths of heavy 
quarkonia, many NRQCD matrix elements are needed. The specific number depends 
on the order in v of the non-relativistic expansion to which the calculation is per- 
formed and on the power counting. At order mv^ and within a conservative power 
counting, S- and P-wave electromagnetic and hadronic decay widths for bottomonia 
and charmonia below threshold depend on 46 matrix elements [51] . More are needed 
at order mv'^ [59l |60l [61] . Order mf corrections are particularly relevant for P-wave 
quarkonium decays, since they are numerically as large as NLO corrections in as, 
which are known since long time [62] and to which the most recent data are sensi- 
tive [63l [1]. NRQCD matrix elements may be fitted to the experimental decay data 
[Ml [65] or calculated on the lattice [661 EH [68]. The matrix elements of color-singlet 
operators are related at leading order to the Schrodinger wave functions at the origin 
[To] and, hence, may be evaluated by means of potential models ^69j or potentials 
calculated on the lattice [T7]. However, a great part of the matrix elements remain 
poorly known or unknown. 

In the matching coefficients large contributions in the perturbative series coming 
from bubble-chain diagrams may need to be resummed [70] . 
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In lattice NRQCD in [32j, the ratio r(T(25) e+e-)/r(T(lS) ^ e+e-) x 
^T{2S)I -^riis) been calculated. The result on the finest lattice compares well 
with the experimental one. 

The imaginary part of the potential provides the NRQCD decay matrix elements 
in pNRQCD. For excited states, they typically factorize in a part, which is the wave 
function in the origin squared (or its derivatives), and in a part which contains gluon 
tensor-field correlators [551 ES EH [57] . This drastically reduces the number of non- 
perturbative parameters needed; in pNRQCD, these are wave functions at the origin 
and universal gluon tensor-field correlators, which can be calculated on the lattice. 
Another approach may consist in determining the correlators on one set of data (e.g. 
in the charmonium sector) and use them to make predictions for another (e.g. in the 
bottomonium sector). Following this line in [55l [58] , at NLO in but at leading 
order in the velocity expansion, it was predicted rhad(xw(2-P))/rhad(Xfe2(2-P)) ~ 4.0 
and rhad(Xfei(2-P))/ rhad(Xfe2(2-P)) ~ 0.50. Both determinations turned out to be 
consistent, within large errors, with the CLEO III data [IJ. One should notice that 



at some order of the expansion in f, the scale ymAqcD start also to contribute in 
pNRQCD jeopardizing in some cases the effective reduction of the nonperturbative 
operators [56] . 

For the lowest resonances, inclusive decay widths are given in weakly coupled 
pNRQCD by a convolution of perturbative corrections and nonlocal nonperturbative 
correlators. The perturbative calculation embodies large contributions and requires 
large logs resummation. The ratio of electromagnetic decay widths was calculated 
for the ground state of charmonium and bottomonium at NNLL order in [7^. In 
particular, they report: V{r]b 77)/r(T(l^) e+e") = 0.502 ± 0.068 ± 0.014, 
which is a very stable result with respect to scale variation. A partial NNLL* order 
analysis of the absolute width of T(IS') e~^e~ can be found in [78j. 




9.2 Electromagnetic transitions 

Allowed magnetic dipole transitions between charmonium and bottomonium ground 
states have been considered in pNRQCD at NNLO in [TH [72]. The results are: 



T{J/^ ^ 7 r/c) = (1.5 ± 1.0) keV and r(T(15) ^^rji,) = (fc^/39 MeV)^ (2.50 ± 0.25) 



eV, where the errors account for uncertainties (which are large in the charmonium 
case) coming from higher-order corrections. The width r{J/ip — > 'yrjc) is consistent 



with [73]. Concerning r(T(lS') ^ ■jrjb), a photon energy = 39 MeV corresponds 




to a r]b mass of 9421 MeV. The pNRQCD calculation features a small quarkonium 
magnetic moment (in agreement with a recent lattice calculation [80]) and the in- 
teresting fact, related to the Poincare invariance of the nonrelativistic EFT [?], that 
Ml transition of the lowest quarkonium states at relative order v"^ are completely 
accessible in perturbation theory [Tlj. 

In the weak-coupling regime, the magnetic-dipole hindered transition T(25') —>■ 
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7^7b at leading order [7T] does not agree with the experimental upper bound [39]. It 
should be still clarified if this is related to the fact that T{2S) system belongs to 
the strong coupling regime or if it is due to large corrections (more relevant in the 
hindered case). 

9.3 Semi- inclusive decays 

The radiative transition T(IS') ^ 7X has been considered in [741175]. The agreement 
with the CLEO data of [76] is very satisfactory when one properly includes the octet 
contribution in pNRQCD [75] . In the same work it is found that the ratios for different 
n of the radiative decay widths r{T{nS) 7^) are better consistent with the data if 
T(IS') is assumed to be a weakly-coupled bound state and T{2S) and T{3S) strongly 
coupled ones [H]. 

In general in exclusive decays and for certain kinematical end points of semi- 
inclusive decays, NRQCD or pNRQCD should be supplemented by coUinear degrees 
of freedom. This can be realized in the framework of Soft CoUinear Effective Theory 
(SCET) [93. 

10 Baryons with two or more heavy quarks 

The SELEX collaboration at Fermilab reported evidence for five resonances that may 
be possibly identified with doubly charmed baryons, see the presentation of Jurgen 
Engelfried at this meeting [2J and [44j • Although these results have not been confirmed 
by other experiments ( FOCUS, BELLE and BABAR) they have triggered a renewed 
theoretical interest in doubly heavy baryon systems. 

Low- lying QQq states may be assumed to realize the hierarchy: m ^ mv 3> Aqcd, 
where mv is the typical inverse distance between the two heavy quarks and Aqcd is the 
typical inverse distance between the centre-of-mass of the two heavy quarks and the 
light quark. At a scale /i such that mv ^ fi ^ Aqcd the effective degrees of freedom 
are QQ states (in color antitriplet and sextet configurations), low-energy gluons and 
light quarks. The most suitable EFT at that scale is a combination of pNRQCD and 
HQET [ini HG]. The hyperfine splittings of the doubly heavy baryon lowest states 
have been calculated at NLO in as and at LO in AQc^/m by relating them to the 
hyperfine splittings of the D and B mesons (this method was first proposed in |47j ) . 
In ESI, the obtained values are: M=* -M= = 120±40 MeV and M=. -M=,, = 34±4 
MeV, which are consistent with the quenched lattice determinations of [HI HQ] and 
the quenched NRQCD lattice determinations of [501 El]- Chiral corrections to the 
doubly heavy baryon masses, strong decay widths and electromagnetic decay widths 
have been considered in [52] . 

Also low-lying QQQ baryons can be studied in a weak coupling framework. Three 
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quark states can combine in four color configurations: a singlet, two octets and a 
decuplet, which lead to a rather rich dynamics [15]. Masses of various QQQ ground 
states have been calculated with a variational method in [53j: since baryons made of 
three heavy quarks have not been discovered so far, it may be important for future 
searches to remark that the baryon masses turn our to be lower than those generally 
obtained in strong coupling analyses. 

For QQQ baryons with a typical distance of the order Aqcd inverse, the form of 
the static, 1/m and spin dependent nonperturbative potentials have been obtained in 
pNRQCD [15]. Up to now only the static potential has been evaluated on the lattice 
[I71I79]. 

11 Gluelump spectrum and exotic states 

Gluelumps are states formed by a gluon and two heavy quarks in a octet configura- 
tion at small interquark distance [99] . The mass of such nonperturbative objects are 
typically measured on the lattice. The tower of hybrids static energies [98] measured 
in lattice NRQCD reduces to the gluelump masses for small interquark distances. In 
pNRQCD [HI [26] the full structure of the gluelump spectrum has been studied, ob- 
taining model independent predictions on the shape, the pattern, the degeneracy and 
the multiplet structure of the hybrid static energies for small QQ distances that well 
match and interpret the existing lattice data. These studies may be important both 
to elucidate the confinement mechanism (the gluelump masses control the behaviour 
of the nonperturbative glue correlators appearing in the spectrum and in the decays) 
and in relation to the exotic states recently observed at the B- factories. The y(4260) 
in the charmonium sector may be identified with an hybrid state inside such picture. 
A complete pNRQCD description of heavy hybrids is still missing. 

12 Production 

Before the advent of NRQCD, colour singlet production and colour singlet fragmen- 
tation underestimated the data on prompt quarkonium production at Fermilab by 
about an order of magnitude indicating that additional fragmentation contributions 
were missing [85]. The missing contribution was precisely the gluon fragmentation 
into colour-octet ^Si charm quark pairs. The probability to form a J/ip particle from 
a pointlike cc pair in a colour octet '^5'i state is given by a NRQCD nonperturbative 
matrix element which is suppressed by with respect to to the leading singlet term 
but is enhanced by two powers of as in the short distance matching coefficient for 
producing colour-octet quark pairs. Introducing the leading colour-octet contribu- 
tions, the data of CDF could be reproduced and this has been an important result 
of NRQCD [85]. NRQCD factorization has proved to be very successful to explain 
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a large variety of quarkonium production processes (for a review see the production 
chapter in [1] and [5]). A formal proof of the NRQCD factorization formula for heavy 
quarkonium production has however not yet been obtained. Recently, there has been 
important work in the direction of an all order proof in |100l llOlj . In particular, it 
has been shown that a necessary condition for factorization to hold at NNLO is that 
the conventional octet NRQCD production matrix elements have to be redefined with 
Wilson lines, acquiring manifest gauge invariance. 

For production, a pNRQCD formulation is not yet existing. In principle to go 
through a further factorization also in production, if at all possible, may reduce the 
number of nonperturbative matrix elements and enhance the predictive power. 

Two outstanding problems exist at the moment in quarkonium production: dou- 
ble charmonium production in e^e~ collisions and charmonium polarization at the 
Tevatron. 

In [102], BELLE reports a{e+e- -> J/ip + Vc) Br(cc > 2 charged) = 25.6±2.8± 
3.4 fb and in [103J . BABAR reports cr(e+e- -> J/ip + rjc) Br (cc ^ > 2 charged) = 
17.6 ± 2.8tli fb- Originally these data were about one order of magnitude above 
theoretical expectations. Recently, with some errors corrected in some of the the- 
oretical determinations, NLO corrections in as calculated in [104j and higher-order 

corrections obtained in |105j . the theoretical prediction has moved closer to the 
experimental one. In [106] . a preliminary estimate of a{e~^e~ J/ip + rjc) including 
the above corrections gives 16.7 ± 4.2 fb. 

Still open is the problem of the BELLE measurement a{e~^e^ J / ip+cc) / a{e^e^ - 
J/ip + X), about 80%, with respect to theory calculation that gives about 10% (see 
[T] for a detailed discussion). 

Charmonium polarization has been measured at the Tevatron by the CDF col- 
laboration at run I with 110 pb"^ [M] and recently at run II with 800 pb"^ [108j . 
The data of the two runs are not consistent with each other in the 7-12 GeV re- 
gion of transverse momentum, pt, and both disagree with the NRQCD expectation 
of an increased polarization with increased pt- For large pt, NRQCD predicts that 
the main mechanism of charmonium production is via color-octet gluon fragmenta- 
tion, the gluon being transversely polarized and most of the gluon polarization being 
transferred to the charmonium. The CDF data do not show any sign of transverse 
polarization at large px- 

A solution to such problem may be obtained in the case in which a nonperturbative 
power counting is valid [551 El]- 

13 Challenges 

For what concerns systems close or above the open flavor threshold, a complete and 
satisfactory understanding of the dynamics has not been achieved so far and a corre- 



15 



sponding general NR EFT has not yet been constructed. Such systems are difficult 
to address also with a lattice calculation. Hence, the study of these systems is on a 
less secure ground than the study of states below threshold. Although in some cases 
one may develop an EFT owing to special dynamical conditions (as for the X(3872) 
interpreted as a loosely bound D*^ + D*^ molecule |109j ). the study of these 
systems largely relies on phenomenological models |110[ 1111] . The major theoretical 
challenge here is to interpret the new states in the charmonium region discovered at 
the B-factories in the last few years. 

Heavy ion experiments use quarkonium suppression as one of the smoking guns 
for quark-gluon plasma formation (cf. e.g. the media chapter in [1]). . To de- 
scribe quarkonium suppression it would be important to formulate an EFT for heavy 
quarkonium in media and to obtain a clear definition of the heavy quark potential at 
finite T . Preliminary studies are ongoing with several approaches |112j . 

With a good control in theory and high statistic data sample available at present 
and future (Super-B factory) experiments, heavy quarkonia may also supply us with 
an alternative way of looking for new physics BSM, cf. [113j and the BSM chapter in 

14 Outlook 

Today NR EFTs and lattice calculations allow us to investigate a wide range of heavy 
quarkonium observables in a controlled and systematic fashion and, therefore, learn 
about one of the most elusive sectors of the Standard Model: low-energy QCD. 

Predictions based on non-relativistic EFTs are conceptually solid, and systemat- 
ically improvable. EFTs have put quarkonium on the solid ground of QCD: quarko- 
nium becomes a privileged window for precision measurements, new physics and con- 
finement mechanism investigations. 

Many new data on heavy-quark bound states are provided in these years by the 
B-factories, CLEO and the Tevatron experiments. Many more will come in the near 
future from BES-III, LHC and later Panda at GSI. They will show new (perhaps 
exotic) states, new production and decay mechanisms and they will be a great arena 
for new EFT tools. 
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